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Geomagnetic Semiannual Variation Is Not Overestimated and Is 
Not an Artifact of Systematic Solar Hemispheric Asymmetry 

L. Svalgaard''^ 



Mursula et al. [2011] (MTLll) suggest that there is a 22- 
year variation in solar wind activity that coupled with the 
variation in heliographic latitude of the Earth during the 
year, gives rise to an apparent semiannual variation of ge- 
omagnetic activity in averages obtained over several solar 
cycles. They conclude that the observed semiannual varia- 
tion is seriously overestimated and is largely an artifact of 
this inferred 22-year variation. We show: (1) that there is no 
systematically alternating annual variation of geomagnetic 
activity or of the solar driver, changing with the polarity of 
the solar polar fields, (2) that the universal time variation 
of geomagnetic activity at all times have the characteristic 
imprint of the equinoctial hypothesis rather than that of the 
axial hypothesis required by the suggestion of MTLll, and 
(3) that the semiannual variation is not an artifact, is not 
overestimated, and does not need revision. 



1. Introduction 

Mursula et al. [2011] noticed that the strongest geomag- 
netic activity during the short 16-year interval 1993-2008 
occurred when the Earth was at southern heliographic lati- 
tudes in 1994 and at northern heliographic latitudes in 2003. 
These two observations are presented as evidence that there 
is an annual variation of solar wind speed at Earth, changing 
phase between the two cycles resulting in a predominance 
of high speed streams from the southern hemisphere during 
the late phase of cycle 22 and from the northern hemisphere 
during the late phase of cycle 23 (although the high geomag- 
netic activity in October 2003 originated from solar activity 
in the southern hemisphere), and that this alternation is 
characteristic of solar cycles in general, providing long-term 
predictability of activity. The phase of the purported effect 
is such that the Earth experiences the fastest solar wind 
when it is north of the solar equator during positive polar- 
ity epochs (magnetic field positive (outward) at the northern 
pole of the Sun) such as during the decline of sunspot cycle 
22 from 1990-1996 and when it is south of the solar equa- 
tor during negative polarity epochs such as the decline of 
sunspot cycle 23 (2000-2008). 

MTLll suggest a novel variation on the axial mechanism, 
namely that this inferred 22-year variation in solar wind ac- 
tivity, coupled with the variation in heliographic latitude of 
the Earth during the year, gives rise to an annual varia- 
tion of geomagnetic activity with opposite phases between 
the declining phases of the cycles, and that these two an- 
nual waves when averaged over several cycles result in an 
apparent semiannual variation. They conclude that the ob- 
served semiannual variation is seriously overestimated and 
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is largely an artifact of a 22-year variation. We examine this 
claim using the full geomagnetic record extending back well 
into the 19th century and find it unjustified. 

2. The Semiannual and Universal Time 
Variations 

As MTLll point out, the semiannual variation is a prob- 
lem of long standing [Mairan, 1733; Broun, 1848; Bartels, 
1932; Mcintosh, 1959; Oliver et al, 2000]. Geomagnetic ac- 
tivity is a very complicated phenomenon. To reduce the 
complexity we commonly resort to study geomagnetic in- 
dices instead. Many indices exist (Ci, u, Dst, AE, IL, ap, 
aa, am, IHV, etc), e.g. Menvielle et al. [2011]. When Julius 
Bartels designed the geomagnetic activity index Kp [Bartels 
et a/., 1939; Bartels, 1949], he found that the i^^ indices for 
each station had such a strong dependence on local time 
that the very non-uniform longitudinal distribution of the 
Kp stations precluded any investigations of the expected 
Universal Time (UT) of the activity as the angle, ^, be- 
tween the direction to the Sun and the Earth's magnetic 
dipole varies both with day of year and with hour of uni- 
versal time, thus introducing different semiannual variations 
depending on Universal Time. Thus it was necessary to re- 
move the local time variation by applying normalization fac- 
tors for each station for each month. The resulting Kp index 
should then have no UT variation at all. The ap index de- 
rived by linearizing Kp is therefore, by design, not suitable 
for a comprehensive study of the semiannual variation of 
geomagnetic activity. Figure 1 shows how strongly the UT- 
variation modulates the semiannual variations at individual 
stations, using the 'raw' IHV index for each station [Sval- 
gaard and Cliver, 2007]. /iJU measures geomagnetic activity 
for an interval around local midnight where the contribution 
(or interference) from the diurnal Sq variation is minimal. 
Calculating the curves using local K indices (converted to 
amplitudes), rather than /i7U yields identical results, as ex- 
pected, because both indices respond the same to the same 
physical processes. 

Mayaud [1967, 1970] utilized a more uniform distribution 
of geomagnetic observatories in both hemispheres to con- 
struct a true planetary activity index. Km, and its linear 
version am. Following Svalgaard [1977], Oliver et al. [2000] 
plotted the am index as 2D-contours as a function of month 
of year and of universal time, to reveal a characteristic 'hour- 
glass' structure. The probability of injection of energy into 
the ring current [O'Brien and McPherron, 2002] and the 
'raw' IHV index [Svalgaard and Oliver, 2007] display, Fig- 
ure 2, the same characteristic hourglass structure, showing 
that the semiannual/UT variation is not just an artifact of 
the am index. New indices based on the same principles pro- 
duce the same hourglass signature, e.g. Figure 3 of Finch et 
al [2008]. Importantly, they show that for the last 50 years, 
a time when the 22-year solar wind-variability claimed by 
MTLll to be particularly well-defined, geomagnetic indices 
show the clear imprint of the equinoctial hypothesis for the 
semiannual variation of geomagnetic activity, [Bartels, 1932; 
Mcintosh, 1959; Svalgaard, 1977; Oliver et al, 2000], rather 
than the absence of the UT-variation one would expect to 
see if the axial effect invoked by MTLll were the dominant 
cause of the semiannual variation. 
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3. Solar Wind Drivers 

Svalgaard [1977] (see also Crooker et al. [1977] and Mu- 
rayama [1982]) has shown that the am index can with good 
approximation be calculated as am — k q(a, DOY, UT)BV^ 
where /c is a scale factor, ^ is a geometrical factor depending 
on the IMF clock angle (a), on day of year {DOY)^ and 
on Universal time (UT), B is the magnitude of the IMF 
strength, and V is the solar wind speed, all at Earth. To 
separate to first order the geometric effects from the intrin- 
sic variations of the solar wind B and U, we can divide the 
geomagnetic activity index under investigation by the prod- 
uct BV'^. Fi gure 3 shows the resulting hourglass (right- 
hand panel) and the variation of the dipole tilt angle. It is 
clear that the semiannual variation is closely regulated by 
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Figure 1. Semiannual variation of Raw /if U for stations 
at different geographic longitudes; southern hemisphere 
stations: reddish color. Upper curves (left-hand scale) 
for stations with longitude near those of the geomagnetic 
poles. Lower curves (right-hand scale) for stations with 
longitudes ^90° away from the poles. Filled symbols 
show observations (roughly spanning the 20th century 
using predecessor stations as needed), while open sym- 
bols show the variation calculated from the Svalgaard 
Function, 6'(^) = 1.175/(1 + 3cos^(^))^/^ [O'Brien and 
McPherron, 2002; Svalgaard, 1977]. All values are nor- 
malized to the mean for the year and repeated for one 
year in the right-hand half of the Figure. 




Figure 2. Left: Probability of ring current injection 
[O'Brien and McPherron, 2002]. Right: Variation of 
"raw" IHV [Svalgaard and Oliver, 2007] for the ^60 sta- 
tions used to calculate IHV. Warm (red) colors signify 
maxima and cold (purple) colors minima. IHV measures 
geomagnetic activity for an interval around local mid- 
night where the contribution (or interference) from the 
regular diurnal Sq variation is minimal. 



the dipole tilt towards the solar wind, regardless of B and 
V. In fact, the hourglass is found for both northward and 
southward IMF as shown in the lower panels. 

Plotting am (normalized to its average) against the tilt 
angle, separately for B and V shows. Figure 4, that the pat- 
tern of semiannual/UT variation of activity with the dipole 
tilt is independent of the intrinsic values of these solar wind 
parameters. The semiannual/UT variation is a permanent 
feature of the interaction between the solar wind and the 
magnet osphere, its functional form independent of the di- 
rection and strength of the IMF and of the speed of the 
solar wind, that is: The dominant semiannual-UT variation 
is a modulation of activity generated by the impact of the 
solar wind, as has been known for decades [Mayaud, 1977; 
Svalgaard, 1977]. The notion, put forward in MTLll and 
elsewhere [Karinen and Mursula, 2006; Love and Gannon, 
2009], that the semiannual variation is 'excessive', an arti- 
fact, due to uneven station distribution, incomplete removal 
of the regular diurnal Sq variation, or other deficiencies in 
deriving geomagnetic indices, is contradicted by the regu- 
larity of the combined semiannual-UT variation observed at 
stations of every longitude, for all solar wind conditions, for 
different indices, as well as for related geophysical phenom- 
ena [Mairan, 1733; Baker et al., 1999]. 

4. Annual Variation of Solar Wind Parameters 

Svalgaard and Oliver [2007, 2010] show how both B and 
V can be determined from the geomagnetic record (the IDV 
and /if y indices). Figure 5 shows 7-rotation running aver- 
ages of the product BV^ comparing the value reconstructed 
from geomagnetic activity directly with in situ observations, 
covering the whole of the space age. It is clear that the re- 
constructed values are a good representation of the physical 
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Figure 3. The variation of am reduced for the influence 
of IMF strength and solar wind speed for the interval 
1963-2003 for all times where we have simultaneous am, 
B, and V data, from Svalgaard and Schulz [2004]. The 
lower panel shows that the existence of the 'hourglass' 
structure does not depend on the orientation of the IMF. 



SVALGAARD: SEMIANNUAL VARIATION 



X- 3 



reality when several rotations are averaged (for individual 
rotations the assumption of balanced northward and south- 
ward fields occasionally breaks down). MTLll suggest that 
the phase of the annual variation of solar wind speed and 
geomagnetic activity changes systematically from one solar 
cycle to another, that the annual variation is largest in the 
declining phase of solar cycles, and that annual maxima are 
located in March during positive polarity periods and in 
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Figure 4. Upper: Relative variation of am (i.e. divided 
by the mean value) as a function of dipole tilt angle, 
^, separately for northward (GSM, blue) and southward 
IMF (red) for odd/even years (diamonds/dots) from Sval- 
gaard and Schulz [2004]. Lower: same but for low/high 
solar wind speed. 
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Figure 5. Seven-rotation running means of BV'^ {B m 
nT and V in units of 100 km/s) calculated from the 
/ifU-/L>U indices, [Svalgaard and Cliver, 2007, 2010] blue 
curve, and determined from in situ observed solar wind 
near the Earth (OMNI: http://omniweb.gsfc.nasa.gov 
red curve). Circles mark the high-speed streams in 1994 
and 2003 that seemed so important to MTLll. Only ro- 
tations with better than 50% coverage were included in 
the calculation. Note the 22-year variation [Cliver et a/., 
2004] in geomagnetic activity resulting in calculated val- 
ues that are slightly too high during the rising phase of 
solar cycles 21 and 23. 



September during negative polarity periods. As shown by 
Svalgaard and Cliver [2007] there is a true, but small (5%), 
second order, variation of geomagnetic activity caused by 
the variation of B with distance to the Sun, as expected 
(their Figure A7 and the insert in Figure 6). This small 
variation, attesting to the validity of the reconstructed solar 
wind parameters, is, of course, independent of the polari- 
ties of the solar cycles, as the orbit of the Earth does not 
change with the cycles, and is not in phase with the varia- 
tions claimed by MTLll. In addition, at any time, the solar 
wind can have enhancements of some duration, especially in 
V as were the cases for the high-speed streams of 1994 and 
2003, or in HMF-polarity anomalies as in 1954 and 1996, 
e.g. Cliver et al. [2004]. 

The central thesis of MTLll 's mechanism for the semian- 
nual variation can be stated as follows: The observed solar 
drivers, the excitation parameters B and U, do not vary ran- 
domly with respect to the Earth's orbit and their variation, 
when coupled with the axial effect, is of sufficient amplitude 
to be the dominant cause of the observed semiannual vari- 
ation. We can put this thesis to a direct test using as our 
measures of the solar driver and of geomagnetic activity the 
long-term variation of BV^ (directly observed since 1965 and 
reconstructed from ID V and /if U before that) and of the aa 
index (which is adequate for monthly values). We consider 
four subsets of the series: rising and declining phases of the 
cycle and for each of those, positive and negative polarity 
(of the northern pole), and then calculate the variations of 
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Figure 6. Upper: Variation of the solar wind driver BV^ 
through the year, normalized to mean values. Bluish col- 
ors for North Pole positive (A > 0; dark blue, squares: 
declining phase; light blue, triangles: rising phase), red- 
dish colors for negative {A < 0; red, diamonds: declining; 
pink, circles: rising). Dashed black curve: all data. Insert: 
variation of normalized BV^ (blue symbols: A > 0; red 
symbols: A < 0; dashed line: all data; full purple curve: 
the inverse square of the solar distance in AU). Lower: 
Variation of aa index; same symbols and colors as for the 
upper plot. 
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BV^ and of aa as a function of month of year. The data for 
this calculation is given in the electronic supplement. It is 
plain from the resulting Figure 6 that during the declining 
phase there is no asymmetry in solar wind parameters nor in 
geomagnetic activity between the first half of the year and 
the second half, contrary to the claims of MTLll. For the 
rising phase, there is no asymmetry either. The claims thus 
fail this direct test. 

As average curves at times may conceal actual variation, 
we show in Figure 7 for individual cycles since 1868 that 
there is no systematic annual variation in the aa index, nei- 
ther for the declining phase or the rising phase, nor for either 
polarity. To compare cycles, the values for each section of a 
cycle were normalized by the average for that section. 

5. The 22-year Variation 

Chernosky [1966] discovered a 22-year variation in geo- 
magnetic activity. The effect is explained [Russell, 1975; 
Svalgaard, 1977] by a combination of the Rosenberg- 
Coleman effect [Rosenberg and Coleman, 1969; Wilcox and 
Scherrer, 1972] varying the dominant polarity of the he- 
liospheric magnetic field and the Russell-McPherron effect 
[Russell and McPherron, 1973], favoring enhanced geomag- 
netic activity when the solar and terrestrial magnetic dipoles 
have opposite directions [Oliver et al, 2004]. This is a 
purely geometrical effect, not dependent on any inherent 
hemispheric asymmetries in solar activity. Echer and Sval- 
gaard [2004] showed that the Rosenberg-Coleman effect is 
only prominent from minimum through the rising phase of 
the solar cycle, so the 22-year cycle is confined to such times 
(the conclusion of Echer and Svalgaard [2004] that the 22- 
year cycle must have another explanation is probably too 
simplistic as they did not appreciate that the enhancement 
of geomagnetic activity is also confined to the same phase 
of the cycle as the Rosenberg-Coleman effect), as is also evi- 
dent in Figure 5. On rare occasions such as in 1954 and 1996 
[Cliver et a/., 2004] this effect (not tied to any inherent en- 
hancement of solar activity, e.g. of solar wind speed) can at 
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Figure 7. Variation of the aa index (1868-2011) sep- 
arately for the declining and rising phases of each so- 
lar cycle and for each polarity, positive or negative. The 
thick, black curve shows the average variation, while thin, 
colored curves are for individual solar cycles. 



those times dominate the semiannual variation, thus mak- 
ing the minor, overall contribution to the average semian- 
nual variation reported by Svalgaard et al. [2002] and Cliver 
et al. [2004]. This (real) 22-year cycle is very distinct from 
and has nothing to do with the variation of solar wind speed 
claimed by MTLll. 

6. Conclusion 

MTLll make two claims 1: that northern and south- 
ern solar hemispheres have a 22-year cycle in systematically 
different, and opposite, activity levels resulting in corre- 
sponding variations of the solar drivers of geomagnetic ac- 
tivity, primarily solar wind speed. When the Earth is north 
(south) of the solar equator, a more active northern (south- 
ern) hemisphere would result in a stronger and faster solar 
wind causing enhanced geomagnetic activity, and 2: that 
the well-established semiannual variation is largely an axial 
effect resulting from the above asymmetry, averaging two 
disparate annual variations offset by six months resulting in 
an artificial semiannual variation. 

We show here that there are no oppositely organized an- 
nual variations in the solar driver of geomagnetic activity 
nor in the observed activity in step with the alternating 
cycle polarities, precluding long-term predictive capability. 
Finally, such a purely axial mechanism does not allow for the 
UT variation which is a well-established and integral part of 
the phenomenon. Consequently, the semiannual variation is 
not 'seriously overestimated' and is not in need of revision. 
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